Dipole emission and coherent transport in random media II. 
Density of states vs. index of refraction 
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This is the second of a series of papers devoted to develop a microscopical approach to the dipole 
emission process and its relation to coherent transport in random media. In this Letter, we deduce a 
relation between the transverse decay rate of an emitter in a virtual cavity and the complex refraction 
index of the host medium. We argue on the possibility of a criterion for inhibition/enhancement 
of spontaneous emission in function of the transition frequency and the correlation length of the 
host scatterers. In addition, we study the radiative/non-radiative nature of the net power emission 
through a microscopical analysis of the scattering events involved. This study reveals essential 
discrepancies with previous interpretations. 
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In a previous work [l[ we computed the average local 
density of states (LDOS) within an isotropic and homo- 
geneous medium. It is the aim of this Letter to use those 
equations to deduce a relation between LDOS and coher- 
ent transport parameters. Recent experiments have 
put forward evidences of such a relation. In [l[ , several 
scenarios were addressed for an emitter seated in either a 
virtual or a real cavity. In this Letter we focus on a vir- 
tual cavity scenario. We will consider the host medium 
made of a set of spatially disconnected point scatterers of 
density p and the emitter as a fluorescence particle seated 
on top of a host scatterer. We will refer to that scatterer 
as host particle. The emitter transition frequency and 
the resonance frequency of the host particle are assumed 
far apart from each other. Dipole induction is completely 
cast in the host particle polarizabiliy. A bulk propagator 
G and electrical susceptibility X can be unambiguously 
defined, where over lines denote rank- 2 tensors. The de- 
cay rate T reads [l[ 
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is the electro- 



In the above formulae, ao = 47ra 3 | !L — ^ 
static polarizability of the host scatterers, being e e (w) 

is 



the single-particle dielectric contrast; a — — — k 

the in-vacuum renormalized single-particle polarizabil- 
ity; p is the spontaneous transition amplitude of the 
emitter and u> the transition frequency, with fcp = lo/c. 



The term — 23?{7^°' ) } accounts for the singularity of the 
radiative part of (r) at r = in vacuum (de- 
noted by the script (0)). The denominator in Eq.{T]) 
amounts to self-polarization cycles over the host parti- 
cle. It gives rise to an effective renormalization of \p\ 2 , 
- , which depends weakly on 7 fac- 
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tors. In a first approach, we can distinguish three con- 
tributions to r" ,M -we drop the scripts on T hereafter, 
namely 
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Eq. (UJ) is the decay rate due to absorbtion within the host 
particle. The other two equations correspond to trans- 
verse and longitudinal decay respectively. Propagating 
modes are determined by the singularities of the inte- 
grands in Eqs. (|5I6| ) in the complex plane, that is, poles 
and branch cuts. Poles are modes fe_L,fc|i satisfying the 
dispersion relations in the effective medium, 

k 2 -k 2 [l + X ±(k)}\ k=k± = 0, 1+X||(*)|k=fc„ = 0. (7) 

The refraction index, n, and the extinction mean free 
path, l ex t, parametrize propagation over long distances. 
That is, they correspond to the longest wave length solu- 
tions for kj_. We will disregard the possibility of propa- 
gating longitudinal modes as these do not affect directly 
the relation between T and the index of refraction. As 
a matter of fact, propagation of longitudinal modes re- 
uires some fine tuning on the arrangement of scatterers 
Because of isotropy, x is a function of k 2 . It is only 
for -C 1 that this dependence is weak, £ being the 
typical correlation length between host scatterers. For 
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a virtual cavity, the cavity radius R can be identified 
with £ as well. It is for /co£ ^ 1 that x admits an ex- 
pansion in powers of k 2 /kg and it suffices to truncate 
the series at order k 2 /k 2 . In doing so, Eq.([5|) presents a 
unique pole k± and the whole transverse propagator is 
characterized solely by n and l ex t, n — + x±(k±)}, 

l~2t — 2fc 3{\/f + XxCATl)}- On the other hand, we can 
distinguish two components in r_|_, namely, 
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where the £ dependent terms are approximate as ex- 
plained below. Further on, invoking causality, n(u>) and 



where the scripts P and NP stand for propagating and 
non-propagating respectively. In Eq.© the 'normalized' 
susceptibility function 3t{ XJ p ^ } plays a role analogous to 
that of the field-strength renormalization factor Z in a 
quantum field theory (QFT). In QFT it is the analytical 
structure of G(k) together with Z in the complex plane 
that determines the mass spectrum of quantum parti- 
cles propagating in space-time [H . Correspondingly, it is 
the analytical structure of G±(k) together with 3t{ Xl p ^ > } 
that determines the spectrum of propagating coherent 
photons in a random medium. In view of the above clas- 

and 2r^ p all together 



sification, we will refer to T° bs , 
as non-propagating decay rate, T 1 



= r° 

abs 
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NP. 



while 2T P is the only propagating component, T p = 2T P 
-compare with other interpretations 15j. It is how- 



ever not clear, except for r° bs , which contribution of 
Y NP is associated to absorbtion by host scatterers and 
which part corresponds to incoherent radiation. That 
needs a deeper microscopical analysis. We proceed to 
in function of the extinction coefficient, 



compute r p 
K = l~ xt /(2ko), and n. 
®{x±(k±)} = n 2 - 1 - 



On the one hand, we have that 
k 2 and 3{xJ-(^-l)} = % nK - On 



the other hand, Eq.© yields 3 " ^{ 
everything together we end up with 
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Taking the approximation e w 
Maxwell-Garnett relation (MG) 



and assuming that 



holds, 
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^-^■y/e. The missing local factor with 



, Eq. lflQj) goes as 

respect to the usual virtual cavity formula [!, [H, 17| im 
plies that only part of the decay rate there is transverse 
and radiative. We will go over this point later on. Pro- 



•=-*/« 



(<0 



\h(r) = -&{r-R) 



r~im 




-O— pfc — 6 — 6 — 6 " "' 

(») 

FIG. 1: (a) Feynman's rules -see also Fights), (b) Diagram- 
matic representation of the total decay rate in a dilute system, 
Eqs. (fT5T^3Tl . 

k(lo) are related through integral equations in cj-modes. 
That is, by Kramers-Kronig relationships (KK). There- 
fore, one can in principle control the variations in 2r^ as 
a function of to and so optimize its value once the varia- 
tions of either n(uj) or l ext (uj) are known [8|]. On the other 
hand, in order to get full control over the total sponta- 
neous emission we have to find similar expressions for T\\ . 
In the limit fco£ 1 we can approximate x\\ ~ XX (^l)- 
MG belongs indeed to this case. In such a situation it is 
possible to give a formula for total r(u;) in closed form 
once a(uj) is known. As the medium surrounding the 
emitter is continuous beyond some distance R = £ in 
application of fco£ <§; 1, Eq.© and the leading order in- 
dependent term of Eq.([5]) approximate the formulae for 
a real cavity in .1] with x(u>) = l _ f "i^( UJ ) > yielding the 
second term in Eq. fTTj) and (see also [l6|) 
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The total decay rate is thus a function of u) and £. Op- 
timization of r with respect to w,( leads this way to 
inhibition/enhancement of spontaneous emission. 

On the other hand, it is expected that not only the 
functional form of e e (w) in host scatterers, but also the 
spatial correlations among these, cooperate towards in- 
hibition/enhancement. For £ ~ k^ 1 the above treatment 
fails as, in general, additional poles in Eq.© enter the 
problem and n, l ex t do not suffice to compute V. Full 
knowledge of x(fc) is needed. Nevertheless, by invoking 
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energy conservation, Ward's identities for vector waves 
Q formulated in a suitable manner may yield integral 
expressions - sum rules over fc-modes - relating real and 
imaginary parts of X\\(k) and x±(k). 

Next, we consider a simple case to study microscopi- 
cally the nature of each of the terms contributing to T 
in Eqs.(@][6]). The setup consists of an emitter on top 
of a host particle in a dilute medium made of spheri- 
cal inclusions of polarizability a. A similar model was 
studied in 1^, In this simple model, the only rele- 
vant correlation is that given by the two-point function 
h(r) = — 0(r — £), £ being the radius of the exclusion vol- 
ume surrounding each scatterer. Therefore, the electrical 

(2) 

susceptibility reads X\\,± — P a + X\\ j_) where 



The first two terms on the right of Eq. (fT6|) equal the 
usual Lorentz-Lorenz (LL) and Onsager-Bottcher (OB) 
formulae in absence of self-polarization 0,0,0], ^ll = 

r o(^) V,r OB = To^) 2 ^ c ryi + |K{ap», 

with e = 1 + Sftjap}. A detailed analysis shows that a 
contribution ^-?it{ap} comes from r" while a contribu- 
tion ^-?R.{ap} comes from 2r~ L . We show this as fol- 
lows. The vertex provided by Eqs. (|13I14|) . which gives 
rise to coupling of longitudinal and transverse modes - 
see FigJ2j6), yields a term common to 2T 1 - and Tu, 
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with cos# = k ■ k' . In the above equations, G±\ are 
the transverse and longitudinal components of the prop- 



A(fc) 
—k 2 ' 



agator of the electric field in free space, G ± (k) 
G\^ (k) — ^§r-; k being a unitary vector along the propa- 

II K Q 

gation direction and A(fc) = I—k®k being the projective 
tensor orthogonal to the propagation direction. In spa- 
tial space, G^{r) = ^1+ [^V® V] At leading 

order in pa, the problem becomes one of single scatter- 
ing. Therefore, the only correlation that matters is that 
of the emitter with the host scatterers, h(r). 

Pictorially, upon self-polarization cycles, Y is given by 
the first three diagrams of FigQJfr). It reads, 
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where To = is the total decay rate in vacuum 

disregarding self-polarization effects. In the limit kgt; <C 
1, we expand T in powers of fco£, 
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In the limit fco£ <C 1, the above expression reduces to 
^■3?{ap}. It equals the contribution of each local field 
factor, ^i 2 - in T^l and -^rf m ^ob- On the other hand, 
2T 1 - includes the term 



d^fc 



l[G^(k)] 2 }, (19) 



which amounts to ^$l{ap}. Its corresponding counter- 
part in Tll and Tob is in the bulk factor s/e. In some 
works only the bulk factor is associated to propagating 
deca y ra te @, [l5| or even to the total transverse decay 
rate [141 ] . Our microscopic approach reveals clear discrep- 
ancies with those interpretations. 

Moreover, one of the common terms given in Eq. (|18|) . 
^3?{ap}, corresponds to T\\ and is associated to 
Ss{x\\/(P a )} i n Eq.flBJ). It is therefore part of the inco- 
herent emission which is scattered away as bare longitu- 
dinal modes couple to bare transverse-radiation modes at 
the cavity surface. The existence of the common term in 
Eq. (|18p is a consequence of reciprocity. However, while 
it gives rise to propagating emission when the coupling 
reads Gii — G±, it gives rise to extinction as read in op- 
posite direction as in Fig[2K&). This longitudinal term 
corresponds to the remaining local field factor in Tll 
and Fob at leading order. Some works @, [TJ, [ijj make 
the identification Vll.ob = 2r^ L OB , which is erroneous 
according to our approach. 

In the presence of absorbtion, the electrostatic polariz- 
ability is complex and additional non-propagating (NP) 
terms in both 2T 1 - and T" contribute to the total decay 
rate. These are, 
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At leading order in fco£, the dominant term is that of 
Eq. (j2Tj) due to absorbtion in the Coulomb-Coulomb cou- 
pling, T NP » Tf s - ~ T ^{aop}^ 

The remaining term r of Eq. ifTTi) contains equal 

contributions from 2r^ p and T|| in agreement with 

Eqs.CEHH. _ 

In absence of absorbtion, the terms in Eqs. (|2QI2ip 

still contribute as considering radiative corrections in a. 
These corrections affect both the host medium and the 
host particle. In the former case, radiative corrections 
show up as a consequence of the coupling between trans- 
verse and longitudinal modes in single scattering events 
[lj]. This is due to the 'internal' correlation function 
h a (r) = 0(r — a) -see FigEfef). Fig^e?) illustrates how, 
for an incoming and an outgoing photon being longitu- 
dinal, the coupling to a transverse internal mode gives 
rise to radiation. Hence, the effective vertex in FigJ^c) 
with a including radiative corrections. This kind of terms 
amount to scattered incoherent radiation. Eq. ([20|) is in 
fact part of Eq.© and so of Eq. (fTT|) . 

Finally, we consider the back-reaction of the host 
medium on the host particle, with electrostatic polar- 
izability cv e -eventually equal to ciq. To this aim and at 
first order, we have to compute the self-polarization dia- 
grams, i.e. 4 th to 6 th diagrams in FigQJfr). For the sake 
of simplicity, we disregard absorbtion and radiative cor- 
rections in a and keep leading order terms in ao -denoted 
by the script (1). We can approximate 
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so that, in absence of self-polairzation, we obtain just 
1 1 1 - T^ B = r (l + l${a a p}) as given above. Plugging 
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Ea.([2"2f into Eq.Q, we get 
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where V£ is the volume of the cavity. Making the identi- 
fication e = 1 + aop, one can verify the agreement of the 
above expression with the OB formula at leading order 
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In summary, we have derived equations which relate 
the spontaneous decay rate of an emitter within a virtual 
cavity, both propagating -Eq. (| 1 0[) — and non-propagating 
-Eqs. (jll|12p . to the mean free path l ex t and the refrac- 
tion index n of the host medium. The last equation being 
strictly valid only in the range of validity of MG. Because 
n and l ex t depend ultimately on uj and £, it is argued 
that variations on T towards inhibition/enhancement can 
be achieved by applying KK relations and Ward's iden- 
tities. Our derivation is based on the approach of []]]. 
That is shown to include both LL and OB at leading or- 
der in p. Microscopical analysis of the scattering events 
involved in the formulae of |l| reveals that local field 



factors present longitudinal decay rate contributions. A 
non-propagating decay rate has been identified. It in- 
cludes both absorbtion and incoherent scattered radia- 
tion. The latter involves coupling of bare transverse to 
bare longitudinal modes which takes place at different 
levels. Those are, at the emitter cavity, at single scatter- 
ers, at the correlation boundary between host scatterers 
and at the boundary of the sample. The latter was not 
considered in this paper. 
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FIG. 2: (a) Feynman's rules, (b) Common vertex derived from 
Eqs. (|13ll4[l which couples longitudinal to transverse modes, 
with \kxk'\ = sin 9. (c) Longitudinal- longitudinal vertex from 
Eq. (|14|l . with \k ■ k \ — cos9. (d) Diagrammatic representa- 
tion of the origin of the coupling of single scattering radiative 
corrections to longitudinal external modes. Integration of in- 
ternal radiative corrections yields the vertex in Fig.(e). 
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